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Applications of perturbative QCD to deeply virtual Compton scattering process require a gener¬ 
alization of usual parton distributions for the case when long-distance information is accumulated 
in nonforward matrix elements of quark and gluon operators. We discuss two types of functions 
parametrizing such matrix elements: double distributions F{x,y;t) and nonforward distribution 
functions P(;{X;t) and also their relation to usual parton densities fa{x). 


PROTON SPIN AND DVCS 


Recently, X. Ji |Q| suggested to use deeply virtual Compton scattering (DVCS) to get information about the total 
quark and gluon contributions to the spin of the proton. Note, that the angular momentum operator can be 
written in terms of the symmetric energy-momentum tensor T°‘0- 


= J d^x , ( 1 ) 

which in QCD can be represented as a sum of quark and gluon parts 

. ( 2 ) 

Hence, the quark Jq and gluon Jg contributions to the proton spin can be obtained from the form factors Aq^g(t), Bq^g(f) 
of the energy-momentum tensor at zero momentum transfer t = 0: 


The functions Aq^g{f), Bq^g{t) are 


'^1,9 ~ 2 [^ 9.9 (*^) + > Jq Jg — 2 ' 

defined by 0 


{p'\Tl^fg\p) = u{p') + Bq,g{t)P^>^a''^-^ 

+ DqJt)Mg>^'^] u{p) 


( 3 ) 


( 4 ) 


where r = p' —p,t = r^ . To get Jq,g, one should know both the non-spin-flip quantities Aq^g(0) related to total hadron 
momentum carried by quarks or gluons and spin-flip amplitudes invisible in deep-inelastic cross sections corresponding 
to exactly forward r = 0 virtual Compton amplitude. However, information about Bq^glt) and Aq^gif), is contained 
in the DVCS amplitude j|]. Anyway, whether the extraction of Bq^g{0) is feasible or not, the studies of DVCS, an 
elastic process exhibiting, in the Bjorken limit, a scaling behavior similar to that of DIS, may be interesting on its 
own grounds 1^,^ (earlier discussions of nonforward Compton-like amplitudes 'y*p —> 'j*p' with a virtual photon or 
in the final state were given in refs. Hi)- 
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PARTON PICTURE: DOUBLE DISTRIBUTIONS 


The kinematics of the process 7 *p ^ ^p' is specified by the initial nucleon momentum p, virtual photon momentum 
q and momentum transfer r = p—p'. We are interested in the situation, when the invariant momentum transfer t = r'^ 
and p^ = rrip are much smaller than the virtuality = —Q^ of the initial photon and the energy invariant p ■ q, with 
the Bjorken ratio I2{p ■ q) = xbj fixed. It is helpful to consider first a (formal) limit = 0 and t = 0 and take 
p as a basic light-cone momentum in the Sudakov decomposition. Since the final photon momentum q' is light-like 
= 0, it is natural to use q' as another basic light-cone 4-vector. Then q = q' — XbjP- Furthermore, in this limit, 
the requirement p'^ = {p + rY = p^ reduces to the condition p • r = 0 which can be satished only if the two lightlike 
momenta p and r are proportional to each other: r = (^p, where C coincides with XBj which satisfies 0 < Xbj < 1- 
The leading contribution in the large-Q^, fixed-a:_Bj, small-t limit is given by DIS-type handbag diagrams 
in which the long-distance dynamics is described by matrix elements (p — r \ 'ipa(0)'y/_i£J(0, z; A)'i/ja(^) j p) and 
(p — r l'tpa(0)"f5j/_tE(0, z; A)'ipa(z) jp), where E(0,z;A) is the usual P-exponential of the gluonic A-field along the 
straight line connecting 0 and z. Though the momenta p and r are proportional to each other r = (p, to construct an 
adequate QCD parton picture, one should make a clear distinction between them. The basic reason is that p and r 
specify the momentum flow in two different channels. For r = 0, the net momentum flows only in the s-channel and 
the total momentum entering into the composite operator vertex is zero. In this case, the matrix element coincides 
with the standard distribution function. The partons entering the composite vertex then carry the fractions Xi of the 
initial proton momentum (—1 < Xi < 1). When x is negative, the parton is interpreted as belonging to the final state 
and Xi is redefined to secure that the integral always runs over the segment 0 < cc < 1. In this parton picture, the 
spectators take the remaining momentum (1 — x)p. On the other hand, if the total momentum flowing through the 
composite vertex is r, the matrix element has the structure of the distribution amplitude in which the momentum 
r splits into the fractions yr and (1 — y)r = yr carried by the quark helds attached to that vertex. In a combined 
situation, when both p and r are nonzero, the initial quark has the momentum xp + yr, while the final one carries 
the momentum xp — yr. In more formal terms, this corresponds to the following parameterization of the light-cone 
matrix elements 

(p - r I V;a(0)zP(0, z; A)V'a( 2 ) |p)U 2 =o = ^(p - r)zu(p) / f 9{x + y<l) 

Jo Jo 

dydx, (5) 

{p-r\tpa{0)l5zE{0,z;A)tjja{z)\p)\:,2^i) = u{p-r)-f5Zu{p) [ [ e{x + y < 1 ) 

Jo Jo 

dydx, (6) 

where z = Though we arrived at the matrix elements (||), (^) in the context of the scaling limit of the DVCS 

amplitude, they accumulate a process-independent information. The coefficient of proportionality between (pz) and 
(rz) is just a parameter characterizing “skewedness” of the matrix elements. The fact that, in our case, ^ coincides 
with the Bjorken variable is specific for the DVCS amplitude. An important feature implied by the representation 
(i,® is the absence of the ^-dependence in the double distributions Ea{x,yA) and Ga{x,y,t). This property and 
spectral constraints ^>0, p>0, x-|-p<l hold for any Feynman diagram Q]. As a result, both the initial active 
quark and the spectators carry positive fractions of the light-cone “plus” momentum p: x + C,y for the active quark and 
(1 — a: — p) -|- (1 — C)y for the spectators. However, the fraction of the initial momentum p carried by the “returning” 
quark is given hy x — yC, and it may take both positive and negative values. 

Taking the limit r = 0 gives the matrix element defining the parton distribution functions fa,a{x), Afa,a{x). This 
observation results in the following reduction formulas for the double distributions E{x,y,t),G{x,y;t): 


pl—X p1 — X 

/ Fa{x,y;t = 0)dy = fa{x) , / Ga{x,y;t ^ 0) dy = Afa{x). 

Jo Jo 


( 7 ) 


The parameterization for the matrix elements given above, results in a parton representation for the handbag 
contributions to the DVCS amplitude: 


|h(p')g'«(p)r)?(c) 

+ ^u{p'){q'f - fq’)u{p)T^iC) + {« ^ a}| 


( 8 ) 


2 





+ie 


)759'«(p) Ta(C) + ^^Hp'hMPmO + {a ^ a} 


where T“(C,t) are the functions depending on the scaling variable 

1 1 


T^iC,t) = - f'dx f " (- 

Jo Jo 


^x-C,y + ie x + C,y 
etc. The terms containing l/{x — C,y ie) generate the imaginary part: 


Fa{x,y-t)dy, 


1 _ 
-lmTp{(,t)= Fa{yC,y;t) 
'x Jo 


dy, 


(9) 


( 10 ) 


with a similar expression for ImT^(^,t). The relation between ImT(^,t) and the double distributions Fa{x,y;t) is 
not as direct as in the case of forward virtual Compton amplitude, the imaginary part of which is just given by parton 
densities /a(C)- Note, that the y-integral in eq.(p^ is different from that in the reduction formula (0), *.e., though 


^a{C,t)= [ Fa{yC,y;t)dy 
Jo 


( 11 ) 


is a function of the Bjorken variable ((, it does not coincide with /o(C)- To get the real part of the \/{x — Cj/ + *e) 
terms, one should use the principal value prescription, be., ReT(C, t) is related to Fa{x, y\t) through two integrations. 


NONFORWARD DISTRIBUTIONS 


Since (rz) = (,{pz), the variable y appears in eqs. ®,(i only in the x + y(^ = X combination, where X can be 
treated as the total fraction of the initial hadron momentum p carried by the active quark. Since C ^ 1 ^od x + y < 1, 
the variable X satisfies a natural constraint 0 < X < 1. Integrating the double distribution F{X — y(, y) over y gives 
the nonforward parton distribution [0j0| 

^min{X/CA/C} 

F^{X-t)= / F{X-yC,y;t)dy (12) 

Jo 


where ^ = 1 — C. The basic distinction between the double distribution F{x,y;t) and the nonforward distribution 
F(^{X; t) is that the latter explicitly depends on the skewedness parameter one deals now with a family of nonforward 
distributions Fi^{X,t) whose shape changes when f is changed. 

The fraction X — f = X' of the original hadron momentum p carried by the “returning” parton differs from X by 
C X — X' = f. Since X changes from 0 to 1 and f ^ 0,1, the fraction X' can be either positive or negative, be., the 
asymmetric distribution function has two components corresponding to the regions 1 > X > ^ and 0 < X < ((. In 
the region X > f, the function X^(X) can be treated as a generalization of the usual distribution function f{x) for 
the asymmetric case when the final hadron momentum p' differs by fp from the initial momentum p. In the region 
X < f the “returning” parton has a negative fraction (X — C) of the light-cone momentum p. Hence, it is more 
appropriate to treat it as a parton going out of the hadron and propagating along with the original parton. Writing 
X as X = X(C, we see that both partons carry now positive fractions Yfp = Yr and yr=(I — Xjrof the momentum 
transfer r. Thus, the nonforward distribution in the region X = Yf f looks like a distribution amplitude '^i^{Y;t) 
for a qq-state with the total momentum r = fp: 

vl/c(r;t)=/ F{{Y-y)f,y,t)dy. (13) 

In terms of the i^-part of the virtual Compton amplitude is 


mo 


X - C + *e 


I 

X-ie 


(X^“(X;<)+X^“(X;<)) dX. 


(14) 


It can be shown that X“’“(X) linearly vanish as X ^ 0 |0]. As a result, the imaginary part is generated by the 
1/(X — ^ -I- je) singularity 
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(15) 


TT >5 

Hence, the integral ^(C?^) in Eq. (0) is equal to i.e., to the nonforward distribution !FQ{X\t) taken at the 

point X = (I- The parameter ^ is present in twice: first as the parameter specifying the skewedness of 

the matrix element and then as the momentum fraction at which the imaginary part appears. As one may expect, 
it appears for X = Xsj = Cj just like in the forward case. Note, however, that the momentum {X — Qp of the 
“returning” parton vanishes when X = (: the imaginary part appears in a highly asymmetric configuration in which 
the fraction of the original hadron momentum carried by the second parton vanishes. Hence, .?y;(C) iu general differs 
from the function /(C). Experimentally, the imaginary part of the DVCS amplitude can be extracted by measuring 
the single-spin asymmetry |^. 
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